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1. Introduction 


One of the most satisfying aspects of the inflationary paradigm is that the enormous structures 
in the present universe actually arose out of minuscule quantum fluctuations in the far past. 
As the spacetime geometry of the universe during inflation was, to excellent approximation, 
that of de Sitter space, this means that the large-scale structure visible today is intimately 
tied to the expectation value, (</> 1 2 ), of quantum fields in de Sitter space. Moreover, through 
observations of the cosmic microwave background and through galaxy surveys, we are able 
to indirectly probe the initial state of the universe. 

Traditionally, the initial state in which the expectation values are calculated is assumed 
to be the Bunch-Davies state, | BD). This state is chosen, firstly, because the corresponding 
Green’s function resembles that of the Poincare-invariant vacuum at short distances and, 
secondly, because (for a massive scalar field) the Bunch-Davies state is invariant under the 
full de Sitter isometry group. Furthermore, the expectation value of the inflaton in the 
Bunch-Davies state, (BD\cj) 2 \BD), leads to the familiar, and observationally corroborated, 
scale-invariant primordial power spectrum of scalar density perturbations. 

However, the assumption of the Bunch-Davies state as the initial state has recently been 
questioned § I- S 0- 0- 10 § I- Several plausible alternatives for the initial state have 
been suggested based on the anticipated behavior of quantum fields at Planckian energies 
© 0 ,0. m 0 0 0 0, [l8| . Because the initial state is modified, these alternatives 
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lead to slightly different predictions for the primordial power spectrum. This is an exciting 
possibility because the signatures of the modified initial state may be within the realm of 
observational investigation p9|, 


21 


23, 


|25| . Nevertheless, one drawback so far has 
been that the assumed behavior of the fields at Planckian energies remains model-dependent. 

In this paper, we take a different approach. We use the fact that energy must be con¬ 
served. Energy conservation does not break down in string theory and it would be surprising 
if we could not assume it during inflation. This suggests a modification of the Bunch-Davies 
state, as we now argue. It is well known that an inertial observer in de Sitter space is immersed 
in a bath of de Sitter radiation emanating from the de Sitter horizon |2(|. The spectrum of 
this radiation (not to be confused with the primordial power spectrum) is, to first approxi¬ 
mation, thermal. Indeed, it can be shown that, in the Bunch-Davies state, the spectrum is 
precisely the thermal Planck distribution characterized by the de Sitter temperature. 

In fact, however, the true spectrum cannot be precisely thermal. Instead, as is known 
from parallel work on black hole radiation, implementing energy conservation modifies the 


spectrum in a definite way [27, 281. One can regard the de Sitter radiation, whose origin is 
quantum-mechanical, as emerging out of virtual pairs of particles, one member of which has 
materialized by tunneling across the horizon. The probability of emission of a quantum of 
energy u) k is then roughly 

T k « exp {-(3uj k ) . (1.1) 

That is, it is roughly thermal. However, it is not exactly thermal because the horizon suffers 
a back-reaction when a particle is emitted, as a consequence of energy conservation. In 


a derivation in which energy conservation is taken into account [29], the effect is that the 
thermal Boltzmann factor is replaced by exp (AS 1 ): 


exp (~/3uJk) -> exp(AS'(a; fc )) , 


( 1 . 2 ) 


where S is the Bekenstein-Hawking entropy of the de Sitter horizon, and A S(uj k ) is the change 
in the entropy when a quantum of radiation with energy u k is emitted. Since the Bunch- 
Davies state resulted in a perfectly thermal spectrum, we may regard this new modified 
spectrum as a consequence of a small modification of the initial state: 

| BD) | BD') . (1.3) 


In turn, we are led to a robust and explicitly calculable modification to the primordial power 
spectrum of inflationary perturbations: 

^P(fc) = J d 3 x e~ i%3 (HD'|0(z)0(O)|HD') • (1.4) 

The modification to the primordial power spectrum thus arises unavoidably by energy con¬ 
servation via a correction to the Bunch-Davies state. 

This paper proceeds as follows. In section 2, we review the different vacuum choices in 
de Sitter space. We also substantiate the following interesting suggestion: the Bunch-Davies 
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state actually appears empty, rather than thermal, to a hypothetical lightlike observer. In 
section 3, we state the results of the tunneling derivation of de Sitter radiation. In section 
4, we present a formalism for obtaining the modification to the Bunch-Davies state from 
the modification to the thermal spectrum. In section 5, we compute the primordial power 
spectrum in this modified initial state. The result, a correction of order ( H/M p ) 2 , while small, 
is universal; as such, it is a model-independent signature of quantum gravity in the sky. 

2. Vacuum States in de Sitter Space 

A well known feature of quantum field theory, and one that makes itself particularly manifest 
in curved spacetime backgrounds, is that the very definition of a particle excitation and hence 
of a thermal state depends sensitively on the choice of vacuum state. In a general curved 
spacetime, there is no canonical or even preferred vacuum state. However, if a spacetime 
admits isometries, and, in particular, a timelike Killing vector field, then this provides a 
natural means of partitioning modes into positive and negative frequency categories and then, 
in line with the standard procedure in Minkowski space, associating these with annihilation 
and creation operators. A vacuum state can then be defined by imposing the condition that 
the state be annihilated by all the annihilation operators. 

In the absence of an everywhere timelike Killing vector, there is no natural choice of 
vacuum state. In that case, one can apply different criteria to motivate particular choices 
of vacuum. A simple choice, where possible, might be to restrict one’s interest to a region 
of spacetime for which a timelike Killing vector does exist, and to use the corresponding 
vacuum state. Similarly, if the spacetime admits an asymptotically Minkowskian region, 
another possibility is to use the natural Poincare vacuum in that region. Alternatively, one 
could demand that the vacuum be annihilated by the generators of some symmetry group. A 
vacuum state can also be deemed unphysical if it fails to satisfy certain criteria. For example, 
if the expectation value of the stress tensor diverges at a nonsingular point in spacetime, such 
as at a horizon, one would consider this grounds to reject the underlying vacuum state. In 
our case, we will require a certain universal form for the particle spectrum detected by an 
inertial observer, a form dictated by the conservation of energy. 

In this section, we review some of the important vacua of de Sitter space, as well as some 
affiliated coordinate systems. 

The Bunch-Davies state 

The Bunch-Davies vacuum state, | BD), is usually defined using planar coordinates, 

ds 2 = —dt 2 + e 2Ht dx 2 = l 9 (—drj 2 + dp 2 + p 2 dfl 2 ) , (2.1) 

(Hr]) z ' ' 

where r/ = —e~ Ht /H is conformal time. See Figure 1. The mode solutions to the massless 
scalar field wave equation, 

u k (r l ,x) = N k (l + ikr 1 )e- ik ^ s , 
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( 2 . 2 ) 



are termed positive frequency modes because they satisfy 

d 


dr] 


u k(jh x) = —ikuk(rj, x) 


(2.3) 


in the infinite past i.e. in the 7] —> —oo limit. An oft-cited motivation for this boundary 
condition is that, in this limit, the physical wavelength of any given mode is arbitrarily 
short compared to the Hubble length, and hence any distinction between de Sitter space and 
Minkowski space should be suppressed. The boundary condition (^d|) is recognized as the 
standard Minkowski space boundary condition, written in de Sitter space. 

Another distinguishing property of the collec¬ 


tion of modes (^2) is that, under an arbitrary 50(4,1) 
de Sitter transformation, the positive frequency modes 
each mix amongst themselves, the negative frequency 
modes mix amongst themselves, but these two classes 
of modes do not mix among each other. This im¬ 
plies, in particular, that the Bunch-Davies vacuum 
state is de Sitter invariant 1 . 

Even so, it is well known that the Bunch-Davies 
vacuum is not uniquely de Sitter invariant; for mas¬ 
sive scalar fields, there is a one-parameter family of 
de Sitter-invariant vacua called a-vacua. These are 
related to the Bunch-Davies vacuum by choosing the 
positive frequency modes to be 

u k(v, x) = A u k (r], x) + B u*_ k (r], x) , (2.4) 

l 



Figure 1: Penrose diagram for de Sit¬ 
ter space. Planar coordinates cover the 
shaded region. The curves are sections of 
constant 77 and constant p. 


with A, B parameterized according to A = y (1 _ eQ+a *^ 
and B = e a A. Here, a is a complex number labeling the a-vacuum; the Bunch-Davies modes 
correspond to B = 0 i.e. to a. = — 00 . A somewhat more transparent expression of the a 


modes is to introduce the antipodal map (r/, x) —► (— 77 , 2 ?) = (x) in terms of which ( 2 T) can 
be written 

u k (x) = Au k (x) + Bu_ k (x) . (2.5) 

Under an arbitrary de Sitter transformation, the positive and negative frequency don’t 
mix, and hence an a-vacuum state, defined by the condition that it lies in the kernel of all 
annhilation operators, is also de Sitter invariant. 


The static state 

Another important state in de Sitter space is the static vacuum, 15). Consider de Sitter space 
in static coordinates: 

ds 2 = -(1 - H 2 r 2 )dt 2 + (1 - H 2 r 2 )~ l dr 2 + r 2 dn 2 . (2.6) 

See Figure 2. The time coordinate here is the proper time of a static geodesic observer at the 
origin and r = H corresponds to the horizon in de Sitter space. 

1 Even though, formally, no de Sitter-invariant vacuum exists for a massless scalar field |^|, by assuming a 

tiny but nonzero mass, we can ignore this subtlety. 
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Using this time coordinate to define positive and 
negative frequencies leads to the static vacuum, so 
called because the generator of time translations, dt , 
is a Killing vector. The physical significance of the 
static vacuum is that the particular geodesic observer 
resting at r = 0 sees it as empty: bk\S) = 0 for all 
k. Of course there is nothing special about this ob¬ 
server, and, indeed, there exists a corresponding vac¬ 
uum state for each of the infinitely many freely falling 
observers that crisscross de Sitter space. The static 
vacuum is invariant under 50(3) x R, i.e. under rota¬ 
tions and time translations, but not under the entire 
de Sitter isometry group, 50(4,1). What seems to 
disqualify the static vacuum as a consistent state in 
de Sitter space is that the expectation value of the 
energy-momentum tensor blows up as one approaches 
the de Sitter horizon [31]. 



Figure 2: Penrose diagram for de Sit¬ 
ter space. Static coordinates cover the 
shaded region. The curves are sections of 
constant t and constant r. 


Relation between the Bunch-Davies state and the static state 


It can be shown that the Bunch-Davies state appears precisely thermal to the static observer 
with temperature 


T = 


H 

277 


(2.7) 


The easiest way to see this is to consider the Bunch-Davies Green’s function, G(x, x'). One sets 
x(t) and x'(t') to lie along a timelike geodesic with proper time r; the observer moving along 
such a geodesic registers a thermal spectrum if two conditions hold. First, the Green’s function 
needs to be periodic under a shift r —» r + /3 (where (3 = T _1 ), and second, the singularities 
of the two-point Green’s function must allow for contour integrals that are consistent with 
detailed balance. In this way, it can be shown that the Bunch-Davies state appears thermal 
to the static observer, and, indeed, by de Sitter symmetry, to all timelike geodesic observers. 

While the Bunch-Davies vacuum is a thermal state, the same is not true of any of the 
a-vacua for a ^ -oo. For B / 0, the combination in ( |2 . 5| ) introduces additional singularities 
for the two-point Green function in the complex plane, which in turn spoils the thermality 
of the state. The a-vacua deviations from a thermal spectrum can be calculated by properly 
accounting for the extra singularities [32, 33, 341, or by using a Bogolubov approach as 
illustrated in Appendix A. 

Perhaps, then, the best way to characterize the Bunch-Davies vacuum is its being the 
unique de Sitter-invariant state that appears thermal to a freely falling observer. Indeed, this 
is an accurate characterization of the Bunch-Davies state, but what is less widely appreciated 
is that there are an infinity of other states that fail to meet this very same set of criteria in 
the most mild of ways. 
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Here is what we mean. Consider a freely falling observer in de Sitter space, using static 
coordinates. Since dt is a manifest Killing vector in these coordinates, we know that 

u 0 = -(1 - H 2 r 2 )^ ( 2 . 8 ) 

is conserved along geodesics. Here r is the proper time and u^ are the components of the 
velocity four-vector (which should not be confused with mode functions). Let’s call the value 
of this conserved quantity —E. Since u 2 = —1 we also have 

= E 2 - (1 - H 2 r 2 ) . (2.9) 

Solving fl2.8| ) and ( |2.9| ) subject to the boundary conditions t s (r = 0) = 0,r(r = 0) = 0 we 
find 

r(r) = (E 2 — 1) sinh(iLr) 

H 

ts{r) = —: tanh _1 (.E tanh(iLr)) . (2-10) 

H 

When E = 1, we find that r(r) = 0 and t s (r) = r, which parameterizes the worldline of the 
static observer who stays put at r = 0. In contrast, notice that for any E ^ 1, t s is invariant 
under r —> r + iriH while r is invariant under r —► r + . 

Consider now the two-point Green’s function for a scalar field, making use of the static 
vacuum, i.e. the zero-particle state as seen by a given E = 1 observer. Since this two-point 
Green’s function when evaluated along a geodesic is a function of (t(r),r(r)), we see that 
the Green’s function is invariant under r —> r + 2iriH~ 1 . Moreover, since in this analysis 
we have not modified the static patch Green function - we have only evaluated it along a 
particular curve - we can be sure that no new singularities have been introduced (in contrast 
to the case with a-vacua). Thus the periodicity of the Green’s function implies that any freely 
falling observer (except for the sole such observer with E = 1) measures a thermal spectrum 
of particle excitations. 

In other words, the static state as defined by the static (E = 1) observer, appears to be a 
thermal state for all other freely falling observers. The only freely falling observer who does 
not see a thermal spectrum is the static observer at r = 0, for whom this is the vacuum. Of 
course, nothing is special about this observer in our analysis; we can do the same analysis for 
every other freely falling observer. Doing so, we conclude that any freely falling observer can 
define a vacuum state with respect to which his or her “Unruh” particle detector will fail to 
detect any particles. The Unruh detectors of all other freely falling observers, however, will 
measure a thermal spectrum with the same temperature. 

There is an intuitive way to understand this. Recall that, for black holes, the state 
that seems empty to the infalling observer looks thermal to the outside observer. Since the 
derivations are identical, this implies that, in de Sitter space, the state that seems empty to an 
observer who falls through another observer’s horizon will seem thermal to the other observer. 
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It follows by de Sitter symmetry that the static vacuum of any observer looks thermal to all 
other observers. 

From the perspective of thermality, therefore, the distinction between the Bunch-Davies 
vacuum and the static patch vacua just described is simply this: All freely falling observers 
see the Bunch-Davies vacuum as a thermal state. All but one freely falling observer sees 
a given static patch vacuum state as thermal. That one special observer, of course, is the 
observer with respect to which the static vacuum is defined. 

Heuristically, this suggests thinking about the Bunch-Davies vacuum state in the following 
way. All observers moving along timelike geodesics see the Bunch-Davies vacuum state as 
thermal. As above, all but one observer sees a given static patch vacuum state as thermal. If 
this one special observer should actually be moving along a lightlike trajectory, then the odd 
man out - the one observer who does not see the static patch vacuum state as thermal - would 
not be among the timelike observers. This suggests that one might think of the Bunch-Davies 
vacuum as that state which appears empty to an observer moving on a lightlike trajectory, but 
which appears thermal to all observers moving on timelike trajectories. From this perspective, 
the Bunch-Davies vacuum and the static patch vacua are all part of a single class of vacuum 
states, each appearing empty to one observer and thermal to everyone else. 

One clue that this heuristic picture is correct is that the light ray emanating from r = 0 
in the limit r/ —> —oo actually travels along the r/ —> —oo time slice. But as i] —> —oo, 
the Bunch-Davies state becomes the ordinary vacuum of Minkowski space. Now, it is known 
that, in Minkowski space, the Poincare-invariant vacuum is the same as the light-cone vacuum 
i.e. the vacuum with respect to the Hamiltonian generating lightlike translations. Thus the 
Bunch-Davies state is the same as the vacuum state for a lightlike observer. To make this 
more explicit, we note that the Bunch-Davies modes become plane waves in the far past. 
Picking an arbitrary spatial direction x , one expands the field operator in the —+ —oo limit 
as 


0 (?y, x, x) = 


r+ °° d 2 k r 1 


dk 7 


\/27t oo 2vr y/2u>k 


{ ( a k x %T e~ iWk1)+ikxX+i%T - ST + h.c.) 

h.c.) } , (2.11) 


+ I a - -iu) k r]-ik x x+ik T -x T i 
' —k x ,k T 


where xt and kr are the transverse position and wave number. The Bunch-Davies state is 
defined by a k £ | BD) = a_ k £ | BD) = 0. But now these modes can be written in light cone 
coordinates x ^ = rj± x: 


(^(x^, x) = 


d 1 | e -i[(uk~kx)x + +(w k +k x )x ]+ik T -x T ^ C ^ 

[(u k +k x )x++(u k -k x )x-]+ik T -x T _|_ ^ c ^ j ^ 12) 


2 tt J — oo 

+ (a 


—k x ,kx 


Since (ui^ ± k x ) > 0, the positive frequency modes defining the Bunch-Davies vacuum also 
correspond to the positive frequency modes defining the light-cone vacuum. Accordingly, 
since one of the light-cone coordinates x^ corresponds to the affine parameter along a lightlike 
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trajectory, the Bunch-Davies vacuum, being equivalent to the light-cone vacuum, also appears 
empty with respect to a fiducial observer traveling along a lightlike trajectory. 

Formally, we can regard the transformation that takes a timelike observer into a lightlike 
observer as an infinite boost. Infinite boosts are not contained in the de Sitter group, which 
is noncompact, but one can include them by extending the de Sitter group to include its 
one-point compactification. Under the de Sitter group, the Bunch-Davies state is special in 
that it is thermal to all geodesic observers. But under the compactified de Sitter group the 
Bunch-Davies state is merely in the same class of states as the static vacua. 


Painleve coordinates 

Painleve coordinates (2^] are a particularly useful coordinate system. The line element is 

ds 2 = — (l — P? 2 r 2 ) dt 2 — 2 Hr dr dt + dr 2 + r 2 dti 2 . (2-13) 


These coordinates resemble static coordinates, but 
have the advantage that the components of the met¬ 
ric and its inverse are finite at the horizon. Moreover, 
constant-time slices correspond to flat Euclidean space. 
In fact, the time slices are equivalent to the constant¬ 
time slices using planar coordinates. See Figure 3. 
Indeed, Painleve coordinates are a hybrid of static 
and planar coordinates, sharing the radial coordinate 
of the former but the time coordinate of the latter: 


2 2 \ 


^Painleve — ^planar — ^static T l n (l Hr 


^Painleve ^static Pplanar ^ 


Ht 


(2.14) 



Physically, Painleve time is the proper time of comov- Fig ure 3: Penrose diagram for de Sit- 

i T-, ,i , ter space. Painleve coordinates cover the 

mg observers, as Painleve time is the same as planar 

shaded region. The curves are sections of 
constant t and constant r. Also shown is 

What is the vacuum associated with Painleve a constant _ time slice split into two. 
time? Because of the off-diagonal metric component, 

the time translation vector dt is different from the normal vector, n, to surfaces of constant 
t. Should one consider modes to be positive frequency if they have positive eigenvalues with 
respect to dt or with respect to n? The answer follows immediately from the Klein-Gordon 
inner product, 


(<£i ,(h) = i f dE a (</>tX<h)- (2-15) 

jt. 

In this expression, dE a = n a dE is the oriented volume element of a spatial slice, £. The inner 
product defines positive norm modes, and, since these are paired with annihilation operators, 
also defines the vacuum state. As the above expression shows, positive frequency should be 
defined with respect to the normal vector, n, not dt . But Painleve time and planar time are the 



same, and hence, so are their time-slices. It follows that the Painleve vacuum is precisely the 
Bunch-Davies state. This is important because the derivation of the approximately thermal 
spectrum in the next section will imply a back-reaction modification to the Painleve state. 
So it must mean that the Bunch-Davies state is modified. 

3. Deviations from Thermality 

In inflationary cosmology, the initial state of the universe is usually taken to be the Bunch- 
Davies state. As explained above, the motivation for this choice of state is two-fold: the state 
is de Sitter invariant and the corresponding Green’s function has a short-distance behavior 
appropriate to Minkowski space. Now, the Bunch-Davies state is a thermal state; an Unruh 
detector for a free-falling observer registers a precisely Planckian spectrum with temperature 
H/2tt. However, there exist quite general reasons why the initial state cannot be strictly 
thermal; the Bunch-Davies state is not physical because it violates energy conservation. To 
see this, consider the closely analogous background of a Schwarzschild black hole. 

At first sight, the spectrum of a black hole’s Hawking radiation is also Planckian. Yet, for 
a variety of reasons, the true spectrum actually departs from pure thermality. First, a purely 
thermal spectrum would have a tail extending out to infinitely high energies. But a black 
hole obviously cannot emit a quantum with more energy than its own mass. Thus there is a 
high-energy cut-off to the spectrum. Moreover, even at energies significantly lower than this 
cut-off, the spectrum must deviate from strict thermality. This is because the “temperature” 
of a black hole is a function of the mass - but the mass is different before and after emission. 

Parallel arguments apply to de Sitter space. De Sitter space has a horizon of radius H ~ 1 . 
This puts an upper bound on the energy of any quantum in de Sitter space: it cannot exceed 
the mass of the “Nariai black hole,” the largest black hole that can fit inside the de Sitter 
horizon. Moreover, as we will see, when the de Sitter horizon emits a quantum of energy, 
it shrinks, much as a black horizon does. The horizon temperature is thus affected by the 
emission of the quantum so that the spectrum has small deviations from thermality even at 
lower energies. This effect become more prominent at higher energies. 

Finally, there is also a more fundamental reason why, in a gravitational setting, one 
does not have a well-defined notion of temperature. Recall that to define a temperature 
for a system, one formally places the system in contact with an infinite reservoir. But, in 
the presence of gravity, no such reservoir exists: any system larger than the Jeans length 
will collapse. One should therefore more properly consider the microcanonical ensemble - 
in which energy is fixed - rather than the canonical ensemble (for which the temperature is 
fixed). 

The above arguments indicate that the Bunch-Davies state, by virtue of being precisely 
thermal, cannot be quite the right state. Rather, the correct state is the one that yields 
the spectrum appropriate to the microcanonical ensemble. Once we have this spectrum, we 
will, in the next section, determine the state that gives rise to it; we will then calculate the 
modification to the power spectrum of primordial density perturbations. 
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Fortunately, there already exists a treatment of de Sitter radiation |29| in a microcanonical 
framework. In this picture, de Sitter radiation arises because a virtual pair is created on the 
other side of the horizon. The positive-energy member of the pair is forbidden outside the 
horizon, because, in static coordinates, the t — t component of the metric has the opposite 
sign there. This positive-energy particle then tunnels across the horizon towards the observer 
where, because of the flip in the sign of the metric, it becomes classically allowed. This picture 
of radiation as a tunneling phenomenon naturally incorporates energy conservation because 
the barrier across which the particle tunnels is determined by the change in the horizon radius. 


One finds [29, 35] that the probability of emission is 


exp(A5) = exp ^—(77 - r\) 


7r 


(3.1) 


where AS is the change in the entropy of the horizon, and r* and 77 are the initial and 


final positions of the horizon. As we will see, at low energies (3T) approximates a thermal 


Boltzmann factor. We note that the form of the expression, as anticipated in [36], is exactly 


the same as for black holes [28, 37]. Another appealing aspect of this result in that context 
is that it is consistent with an underlying unitary quantum theory [38]. 2 


To see that r* differs from 77 , consider a single particle emitted by the de Sitter horizon. 
Before emission, the spacetime is empty and is described by the line element of de Sitter 


space, (273). Thus the initial horizon radius is simply 


r, = H 


-1 


(3.3) 


After emission, the spacetime is not pure de Sitter space because it contains a quantum 
of energy uj. For simplicity, consider emission in the s-wave so that the particle is really a 
spherical shell. The spherically-symmetric spacetime with energy uj is Schwarzschild-de Sitter 
space. This has the line element 

ds 2 = —(1 — H 2 r 2 — 2 Guj/r)dt 2 + (1 — H 2 r 2 — 2Guj/r)~ 1 dr 2 + r 2 dfl 2 . (3.4) 

The final location of the horizon is determined by setting g rr to zero and solving the resulting 
cubic equation for its largest root. We find 

2 (l ( \Ih - (GuHf W 

r/ = 7773 008 (3 aKl “ [ J I • (3 5> 

2 Were there a description of de Sitter radiation in a unitary theory of quantum gravity, one would expect 
the emission rate to be given by the square of the amplitude times the phase space factor. The latter is given 
by summing over the e Sfinal final states and averaging over the e Sinitial initial states. Thus we would expect 

^ g'Sfinal 

T = lamplitudel x (phase space factor) ~ —=-= expfAS 1 ) , (3-2) 

g'-’initial 

in accord with our expression. 
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where the arctan takes values between 7t/2 and 7 r. Equation ([Ll|) then implies an emission 
rate of 


exp 


7T 


H 2 G 


- cos 
3 


1 

— arctan 
3 





-GuH 


- 1 


(3.6) 


This expression can be rendered more tractable by considering the low-energy limit GujH <C 1. 
Then 

27T (_ ljjH \ 


exp 


--U.U + 


87 tM 2 ) j 


(3.7) 


where we have expressed Newton’s constant in terms of the Planck mass: M~' 2 = 87 tG. If 
we neglect the u>H/(8ttM 2 ) correction, we find a probability that is precisely the Boltzmann 
factor for emission at temperature H/2n. Had there been no correction term, elementary 
statistical mechanics would then have implied a spectrum with Planckian occupation numbers, 
precisely the spectrum detected by a free-falling Unruh detector in the Bunch-Davies state. 

The correction term can be regarded as a consequence of self-gravitation or, equivalently, 
of back-reaction or energy conservation |27,|37]. In the next section, we will infer the corrected 
initial state that corresponds to the corrected spectrum. Note that we can already anticipate 
a modification to the primordial power spectrum. Prom ([h?]), we can attribute an effective 
temperature 

H <■ ' (3.8) 

whereas at 


TeS ~ 2 vr V 1 87 tM 2 ) 

so that, at typical energies, u ~ H, we have a modification of order ( H/M p ) 
Planckian energies we find a modification of order H/M p . 


4. Determining the Correction to the Bunch-Davies State 

We will now interpret the corrections to the thermal spectrum in terms of a small modification 
of the initial state. Let us first explain the setup. In canonical language, the probability for 
pair creation of quanta of energy ujk is 


T k = 


(out | b k b-k | in) 


(out | in) 


(4.1) 


where the operators b annihilate the out-vacuum, i.e. fr/, : |out) = 0. We need to determine 
what |in) and |out) refer to in the tunneling calculation. Now, that calculation was performed 
using Painleve coordinates 3 [29, | 


The tunneling probability is the probability to go from 


empty de Sitter space to a spacetime containing a pair of particles, one on either side of the 
horizon, only one of which is detected as de Sitter radiation by the observer at r = 0. This 
suggests that the out-state should be related to the static vacuum | S), as defined by an r = 0 


3 See |27| for earlier work incorporating back-reaction using a Hamiltonian formalism. 
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observer. But, more precisely, because Painleve coordinates cover the full planar patch, the 
out-state actually splits (see Figure 3) into a tensor product: 

|out) = |I) ® |II) , (4.2) 


with the state in the static region I corresponding to the static vacuum, |I) ~ |S). That |I) 
corresponds to the static vacuum can be seen in a couple of ways. First, it is suggested by 
the near-thermal result of the path integral. And second, in region I, the timelike Killing 
vectors <9t static and <9t Painleve agree, by ( 2.14|) 4 . Then since the quanta in the out-state have 
to be created in pairs, and because only one member of the pair ends up on each side of 
the horizon, only one of the creation operators acts on the static vacuum. The other one 
acts on II, on the other side of the horizon. That is, the amplitude is for the transition 
|in) —> b^b^l^\out). Thus the probability of creating two particles in the out-state is the 
same as the probability of detecting one particle in the static vacuum. 

Furthermore, the in-state is, to first approximation, just the Painleve vacuum. We argued 
in section 2 that that was just the Bunch-Davies state, | BD). However, were the in-state 
really the Bunch-Davies state, we would have obtained a precisely thermal spectrum. The 
not-quite thermal spectrum suggests instead that the in-state is a slightly modified Bunch- 
Davies state, | BD'). Indeed, the effective geometry seen by a self-gravitating shell is not pure 
de Sitter space, but rather Schwarzschild-de Sitter space, This has a slightly different 

Painleve metric and is therefore associated with a slightly different state. 

In summary, we have three states: |out), | BD), and \BD'). Moreover, 


(out| b k b-k |BD) 
(out|BD) 

(out| b k b_ k |BD') 2 
(out|BD') 


e~ 0Uk 


gA S(w k ) 


(4-3) 


where the first expectation value is the standard thermal expression, and the second expec¬ 
tation value was derived through the tunneling calculation. Let us now use these expressions 
to determine the precise relationship between the Bunch-Davies state and the modified state 
in terms of Bogolubov coefficients. 


4.1 A general Bogolubov approach 

Consider a real quantum scalar field in the planar patch of de Sitter space. Expand the scalar 
field in a complete set of orthonormal momentum eigenmodes: 

4 >( x ) = <^( x ) i ( 4 - 4 ) 

k 

4 In principle one should be able to determine the structure of the |out)-state unambiguously from the path 
integral derivation; it would be nice to show this formally. 
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with k = \k\. The time coordinate is assumed to be Painleve time ( |2.13| ), which itself is 
equivalent to planar time. Define three different expansions of the same quantum field: 

4> k (x) = a k u k (x) +a}_ k u*_ k (x ) 

= a' k u k (x ) + a_ k u'*_ k (x ) 

= h v k (x) + b ] _ k v*_ k (x) . (4.5) 

Also define three different vacuum states, each corresponding to an “empty” state in terms 
of the associated annihilation operators: 

a k | BD) = 0 , a! k \BD’) = 0 , 6fc|out) = 0 . (4-6) 


We have in mind that these three states are the Bunch-Davies state, the modified Bunch- 


Davies state, and the out-state. As in (4.2), the out-state can be written as a tensor product 
of states in regions I and II; the vacuum state in region I is just the static vacuum. 

The particle content of our three states will generally be nonzero with respect to each 
other; Bogolubov transformations mix up the positive and negative frequency modes. Ex¬ 
pressing the linear decomposition for the mode functions v k (x) in terms of u k (x) or u' k (x) we 
have 


v k (x) = ^2 \a jk Uj(x) + (3 jk u-j(x)* 
j 

= [“ifc u i(x ) + Pik u '-i( x Y ■ 

i 

Similarly, the different creation and annihilation operators are related via 

h = [ & j*> a i ~ P*jk a -j 

3 

E\ ~f* / nl* /1 

y^-ik ® i Pik ® —i 

i 

Proper normalization of the scalar field requires that 



k k 


Pjk Pik Qtjk'j ~ y ] {^ik Pjk Pik ^jk'j ~ 0 * 

k k 


(4.7) 


(4.8) 


(4.9) 


A general Bogolubov transformation, viewed as a matrix, can have nonvanishing off-diagonal 
components. However, for our present situation, because the spatial mode-functions are 
defined on the same spatial slice for all three states 0 , we can choose an orthonormal basis 
such that the Bogolubov coefficients are diagonal in k: 

Ot-kl = OLk^kl i Pkj = Pk^kj i |dfc| \Pk\ 1 

&kl = “fc $kl , Pkj = Pk^kj , | Ol'kV ~ \Pk\ 2 = 1 ' (4-10) 
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We are interested in solving the following problem: Suppose we know the Bogolubov 
coefficients ak, (3k and a' k , (3' k , relating the | BD) and | BD') states to the |out) state. How do 
we then determine the Bogolubov coefficients between the | BD) and | BD') states? In other 
words, we want to express the Bogolubov coefficients 07 and (3k, defined by 

a' k = a* k a k - P k aft (4.11) 


in terms of the coefficients ak, f3k and a! k and / 3' k . Since the final result will be important 


for the rest of the paper, let us be as explicit as possible. Multiplying (f4.8|) , by a' k and its 
complex conjugate by f3 k , and subtracting, we see that 


O' k 


= a' b k + ft b ] 


k —k 


(4.12) 


Using (fO|), we find 

a'k = {ftkft ~ ftPk'j a k + o! k ft + ftWk'j o)_ k ’ 


so that, by (4.11), 


(4.13) 


a* k = a' k a* k -ftp k 

—flk = —ftft + ft^k ■ (4-14) 

As a check, note that when | BD) = \ BD') the Bogolubov coefficients ak and (3k become 
trivial. Once the Bogolubov coefficients a' k , j3' k , ak and (3k are known, these expressions allow 
us to immediately determine the Bogolubov coefficients relating the | BD) the | BD') state. 


4.2 Relating the two in-states 

The Bogolubov transformation guarantees that both the | BD) and the | BD') in-state can be 
realized as squeezed out-states: 


| BD) = S |out) , | BD') = S' |out) 

where S and S' are unitary squeezing operators: 

S = II ( X “ l^fcl 2 ) 3 ex P 


(4.15) 


(4.16) 


The other squeezing operator S' is obtained by replacing 7 k with ~j k . Demanding (Tt) and 
expressing the annihilation operators for | BD) and | BD') in terms of bk and &j), we find 


~ = = _ft 

fk — - ) Tfc — -/ 

a k u' k 


(4.17) 
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This is easiest to show by writing bk as in the Bogolubov transformations. We can also 
show (see for instance |l9[ ^(J) that 

(out | b k b-k JBD) 

(out|BD) 

(outI b k |BD 7 ) 

(out|BD') 

In words: the two 7 / 7 s are the probability amplitudes for pair-creation in the Bunch-Davies 
and modified Bunch-Davies state. 

Now recall that we defined the Bunch-Davies and modified Bunch-Davies state through 
the absolute value of exactly these squared expectation values in fl4.3|) . As a consequence the 
absolute values squared of the Bogolubov ratios 7 *, and 7 )) are 


7 k = 

7 k = 


m 2 = e~^ , |^| 2 = e A5K) - 


(4.19) 


_ ft 


(4.14). We find 


Given these expressions for % and yj) we are now in a position to calculate 7 *, = —using 

1 - 7fc/7fc 


7 k = Ik 


a k 1 -7fc7 fc 


(4.20) 


Now, multiplying one set of operators, say { 6 fc}, by overall phases clearly does not affect 
the choice of state. We can use this rescaling freedom to make a' k real. (Indeed, by further 
rescalings, we can also force af. to be real.) But to proceed, we also need to know something 
about the phase of 7 G In fact, it turns out that neither j' k nor 7 % have any nontrivial phases; 
they are both real and positive. 

There are a couple of ways to see this. Generally, probability amplitudes for tunneling - 
which is what the 7 / 7 s represent - are real. For example, in nonrelativistic quantum mechanics, 
the probability amplitude for a particle to tunnel across a potential barrier V ( x ) is 


i/>(xii) 

i/>(x 7 ) 


exp 



(4.21) 


where the momentum, p, is purely imaginary: p = y / 2m(£’ — V (x)). Thus there is no 
nontrivial phase in the amplitude. 

Moreover, one can show that if 7 *, is real and positive, then 7 I, must be real and positive 
as well. That is, if the amplitude for tunneling without back-reaction has no phase, then 
neither does the amplitude when back-reaction is incorporated. To see this, imagine the 
particle or shell of energy being made up of a large number, N, of smaller noninteracting 
sub-shells, each of energy u^/N. As N becomes infinite, the constituent shells have negligible 
energy and Hawking’s thermal formula, which neglected back-reaction, becomes exact. The 
probability amplitude for emitting one of these sub-shells is then exp(— ^/3bUk/N), where (3 is 
the inverse de Sitter temperature. But now the probability amplitude for the finite shell is 
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just the product of the probability amplitudes for the sub-shells, taking into account that /3 
changes infinitesimally with each emission. Hence 

N 

^ = \\e~^ Uk/N = e-^P Uk/N -f e~*I pdu . (4.22) 

Writing dM = —dcok and invoking the first law of thermodynamics, f3dM = dS, we obtain 


7k = 


(4.23) 


This is perhaps the easiest way to see that the tunneling probability had to take the form 
exp(AS'). But in addition one learns that, since the infinitesimal probability amplitudes had 
no phase, 7 I is real and positive. As a final check, we confirm explicitly that these arguments 
are borne out by considering the particle spectrum of the a-states, for which the result is 
known. This is done in the appendix. We indeed find that 7 % has no nontrivial phase. Then, 
by the above arguments, neither does 7 ^. 


Thus we can write (4.20) as 


7fc = 


7 k 


7 k 


1 - 7fc7 fc 


Plugging in (4. 19|) yields 


g 2 g2 

7fc j x 

l _ g-^/^gi As ( u k) 


(4.24) 


(4.25) 


As a check, we note that 7 *. vanishes in the limit A S(u k ) = —(3u k . i.e. the corrections due to 
back-reaction vanish when —» 0 . 


In (|3.7[) we found that 




(4.26) 


Expanding the exponential in (4.25) to lowest order in u^/Mp, and using f5 = 2n/H, we find 
that 


7 k 


u k H 


(4.27) 


_ 1 8 Ml • 

The energy uo k corresponds to the physical energy of the emitted quantum as measured by a 
static observer at r = 0. In terms of comoving momentum, which is the natural momentum 
label in planar coordinates, the physical energy u k is related to comoving momentum by the 
scale factor: oc it)' We will need this relation in the next section, when we calculate the 

primordial spectrum of inflationary perturbations using the modified Bunch-Davies state. 


5. Primordial Spectrum of Inflationary Perturbations 

We have argued that the standard Bunch-Davies vacuum state needs to be modified because 
it is inconsistent with energy conservation. Rather, the spectrum of perturbations should 


- 16 - 













be calculated instead in the modified state \BD'). Now, the standard primordial power 
spectrum of inflationary perturbations is directly proportional to the massless scalar field 
spectrum; the results we obtain for the scalar field therefore also pertain to the spectrum of 
inflationary density perturbations, as is usual in such calculations. 

Consider then an arbitrary vacuum state |X). Let |X) be annihilated by all the anni¬ 
hilation operators a' k , and call the corresponding Fourier-transformed mode functions v! k {rj). 
The equal-time variance of the scalar field in |X), 

(X\4,(i)m\x) = / (0r Kl 2 , (5.1) 

defines the power spectrum, (X\4> 2 \X) = f din k P(k), of scalar field perturbations: 

2tt 2 r 

-grP(k) = J d 3 xe- ik -* (X\^x)m\X) = K I 2 • (5.2) 

The well-known scale-invariant result arises when the variance is computed in the Bunch- 
Davies state: |X) = | BD). We denote this by Pbd = 777 \uk\ 2 - Now let |X) = \BD'), and 
call P(k ) the power spectrum in the modified state. From the Bogolubov transformation 

u k(v) = a k u k(ri) + Pkutiv) , (5-3) 


it is straightforward to find the relation between the two power spectra. We know that we 
can set ak to be real. Moreover, in the previous section we showed that 7 k = —fd k / a k is real. 
Also, write Uk = \u k \e lS . Then we conclude that 

P(k) = 2 (l - 27 fc cos(2<5) + 7 1) Pbd ■ (5.4) 

^ k 


When calculating the power spectrum of inflationary perturbations one is usually instructed 
to evaluate the spectrum for different k at the time of horizon crossing, corresponding to 
k = a(t)H. The reason this gives the right answer is that, after horizon crossing, the mode- 
function very quickly approaches a constant. Furthermore, the Bunch-Davis mode-functions 
also have the property that the phase 5 is fixed and independent of k, after imposing k = aH. 
We can determine 5 from the expression for the Bunch-Davies modes, (2.2). Also, since 
7 ^ <S 1 , we can write 

P(k) « (1 - 2'y k cos(25))\ k=aH P B d ■ (5.5) 


Now we simply plug in the expression for 7 k , ( |4.27|) , to find the leading correction to the 
primordial power spectrum from picking the modified Bunch-Davies vacuum \BD'). We also 
have to impose the condition k = aH (which leads to cos(2<5) ~ — 1 ). As noted above, this 
condition tells us to evaluate the expression for the power spectrum at the time the physical 
momentum equals the Hubble scale H. Since the expression for 7 ^ is in terms of the physical 
energy scale, we must set = H. Thus we find 


P(k) 


Pbd 


1 + 


Ae 71 



(5.6) 
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This is our final result. We have introduced a ^-dependence because, in a generic slow- 
roll inflationary model, H depends on k. The deviation from unity of the term in brackets 
represents the contribution of the most universal quantum gravity effect - back-reaction due 
to energy conservation - to the primordial power spectrum of inflationary perturbations. 


6. Discussion 


Our calculation shows that the inclusion of back-reaction results in a mild breaking of de 
Sitter invariance with a corresponding mild effect on the CMB power spectrum. This can 
be seen by noting that the Bogolubov coefficients relating the modified Bunch-Davies state 
to the standard de Sitter invariant vacuum depend on the momentum scale cok- Because the 
set of all de Sitter invariant states is paramerized by /c-independent Bogolubov coefficients 
(the a-vacua) this necessarily implies that the modified Bunch-Davies state breaks de Sitter 
invariance. Moreover, as the conservation of energy is mandatory, this observation could 
perhaps have other important consequences, such as for the stability of de Sitter space. 

The correction we computed can intuitively be understood as coming from a slightly mod¬ 
ified effective de Sitter temperature. By instead regarding the effect as a modification of the 
Bunch-Davies vacuum we were able to compute the exact coefficient of the correction. Gener- 
ically, modifications to the initial state during inflation, as studied in [||, 16|, give rise to a 


quite specific modification to the power spectrum: an oscillatory correction with amplitude 
on the order of ( H/M ) where M is the scale at which the modification is implemented. The 
fact that the power spectrum correction we compute here is not of that form, but is instead 
a shift of order (H 2 /M 2 ) - exactly what one expects from including higher order operators 
in a bulk effective field theory |hi|] - suggests that there may be a bulk effective field theory 
description of the modification we’ve studied. In particular, a correction of this magnitude 
generically arises from a Planck-scale-suppressed dimension 6 irrelevant operator; it would 
be interesting to study the exact details and the interpretation in effective field theory of the 
operator responsible for this gravitational correction 5 . 

Our final formula, (|5.6| ), obviously implies a tiny correction, already highly constrained 
by H/Mp < 1CD' 5 . Nevertheless, the result is interesting for a couple of reasons. First, it is 
unavoidable, originating as it does in energy conservation. And second, the correction can be 
calculated explicitly without requiring any knowledge of new Planck scale physics. In that 
sense, it is a universal imprint of quantum gravity in the large-scale structure of the universe. 
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A. Particle Spectrum in the a-States 


Rather than choosing the Bunch-Davies state, which gives rise to a thermal spectrum for a 
free-falling timelike observer, one could just as well have picked one of the de Sitter-invariant 
a-states instead. As an illustration of our general Bogolubov method we shall calculate here 
the spectrum detected by an inertial observer in an a-state. We will reproduce the known 
result [^, [33|, [34|], a result usually obtained by carefully analyzing the singularity structure 
of the Wightman Green’s function. 

Now, by definition, the operators, a'k, that annihilate an a-state are obtained from the 
Bunch-Davies operators, a*,, by a Bogolubov transformation, 

a'k = a* k a k - (5* k a ] _ k , (A.l) 


where the Bogolubov coefficients are ^-independent: 

1 

a k = 7= = = 

s/1 - e a+a 
e a 

^ ~ s/1 - e a+ot * ' 


(A.2) 


Here a is a complex number, with Re(a) < 0, that labels the a-state (and is not to be 
confused with the Bogolubov coefficient a k )\ Re(a) = —oo corresponds to the Bunch-Davies 
state. 

The operators, b k , that annihilate the out-state are themselves related by a Bogolubov 
transformation to the Bunch-Davies creation and annihilation operators: 


b k = a* k a k - (3* k a}_ k . 


(A.3) 


Now, multiplying b k by an overall phase obviously does not affect the choice of vacuum state. 
Use this freedom to fix the phase so that a k is real. Next, express these operators in terms 
of the a-state operators: 


7 / nf* / T 

bk — cv k a k — Pk a —k • 


(A.4) 


Then, by inverting (A.l), 


Clk — Oi k ® k T P k ® —k 


and substituting into (A.2), we conclude that 


(A.5) 


Or, using (|aJ), 


otk = u*k<Xk - P*kPk 
~P'k* = Plat - alPl . 


-/* 

a k 


Pk 


1 

s/l - e a+a * 

1 

s/1 - e a+a * 


r-4 - e a ~P* k 
P*k-e a *a*k 


(A.6) 


(A.7) 
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5/* 

The spectrum is controlled by the absolute value of the parameter fi' k = Now we use 

the fact that, for the thermal Bunch-Davies state, | 7 fc | 2 = e~^ k . Then a little algebra yields 

1 + e°* 


' 12 _ e -/ 3 u k 


\W = e 


% 1 


l + e“7,* 


(A. 8 ) 


In section 4.2, we argued that, on general grounds, 7 *, is purely real. This allows us to write 
7 k = exp (—^/3uk), and we obtain 


ftl 2 = e-* 8 "* 


l _|_ g Q +2^ a 'fc 

1 _)_ e a -2P U k 


(A.9) 


in perfect agreement with previously derived results |33|]. (Had we not known that 7 & was 
purely real, we would have learned that by matching our result with that in the literature.) 
Note incidentally that, for Re(a) 7 ^ — 00 , the spectrum is not thermal; in the limit ^ —> 00 , 
for instance, we find | 7 fc | 2 = e a+a , corresponding to constant occupation numbers. 
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